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We present a numerical method to compute path integrals in effective SU(2) Yang-Mills theories. 
The basic idea is to approximate the Yang-Mills path integral by summing over all gauge field 
configurations, which can be represented as a linear superposition of a small number of localized 
building blocks. With a suitable choice of building blocks many essential features of SU(2) Yang- 
Mills theory can be reproduced, particularly confinement. The analysis of our results leads to the 
conclusion that topological charge as well as extended structures are essential elements of confining 
gauge field configurations. 
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I. INTRODUCTION 

A common approach to get some insight regarding the 
mechanism of confinement is to approximate the Yang- 
Mills path integral by restricting the integration to cer- 
tain subsets of gauge field configurations. By considering 
such effective theories one might find out, which gauge 
field configurations are responsible for confinement and 
which arc not. Examples are ensembles of singular gauge 
instantons (c.f. e.g. [l|), ensembles of calorons (c.f. e.g. 
[3,0]), ensembles of regular gauge instantons and ensem- 
bles of merons 0, [|[ or the removal of center vortices 
in lattice calculations (c.f. e.g. @). Some of these ap- 
proaches have analytical aspects but most of them ex- 
tensively resort to numerical methods. 

In this work we generalize the ideas and techniques pre- 
sented in [J, [f| . We study different classes of gauge field 
configurations, especially their importance with regard 
to confinement, by applyinga numerical method called 
pseudoparticle approach @,H, [|[ . We demonstrate that 
with a suitable choice of building blocks many essential 
features of SU(2) Yang-Mills theory can be reproduced, 
particularly confinement. 

This paper is organized as follows. In section |TT] we ex- 
plain the basic principle of the pseudoparticle approach. 
We also discuss numerical issues and compare the method 
to lattice gauge theory. In section IIIII we show how to 
calculate certain observables in the pseudoparticle ap- 
proach: the static quark antiquark potential at zero and 
at finite temperature, the topological susceptibility and 
the critical temperature of the confinement deconfinc- 
ment phase transition. We present numerical results, 
which are in qualitative agreement with lattice results. In 
section lTVl we apply the pseudoparticle approach with dif- 
ferent types of building blocks, to study different classes 
of gauge field configurations and their effect on confine- 
ment. By doing that we try to determine properties of 
confining gauge field configurations. Our findings indi- 
cate that topological charge as well as extended struc- 
tures are necessary to obtain confinement. In section FVl 
we summarize our results and give a brief outlook regard- 
ing future research. 



II. THE PSEUDOPARTICLE APPROACH 

A. SU(2) Yang-Mills theory 

In this work we consider Euclidean SU(2) Yang-Mills 
theory. The action is given by 
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(g is the dimensionless coupling constant) and field 
strength 

C = d,Al-d„A; + e abc AlAl (3) 
(ei23 = 1)- The topological charge is defined by 

d^xq (4) 
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with topological charge density 
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B. The basic principle of the pseudoparticle 
approach 

The basic idea of the pseudoparticle approach is to re- 
strict the Yang-Mills path integral to those gauge field 
configurations, which can be represented as a linear su- 
perposition of a small number of building blocks, typ- 
ically around 400. The building blocks are gauge field 
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configurations, which are localized in Euclidean space- 
time. We refer to these building blocks as pseudoparti- 
cles, where the term pseudoparticle refers to the fact that 
the corresponding gauge field configurations are also lo- 
calized in time in contrast to ordinary particles, which 
are only localized in space. In our context a pseudoparti- 
cle is not necessarily a solution of the classical Yang-Mills 
equations of motion. 

Every pseudoparticle or building block has certain pa- 
rameters, a position, an amplitude and a color orienta- 
tion, which uniquely define the resulting gauge field con- 
figuration. It is important to stress that in general a 
sum of such pseudoparticles is not even close to a classi- 
cal solution, as it is the case e.g. in instanton gas models 
(c.f. e.g. [H). The pseudoparticle approach is supposed 
to describe full quantum physics and not only certain 
semiclassical corrections. 

The integration over all gauge field configurations in 
the path integral is replaced by an integration over the 
amplitudes and color orientations of the pseudoparticles. 
When considering a spacetime region containing a finite 
number of pseudoparticles, this is an ordinary multidi- 
mensional integral, which can be computed by Monte- 
Carlo methods. 

The starting point of this work has been 0, 0| . How- 
ever, there are two important generalizations: 

• We do not restrict our approach to regular gauge 
instantons and merons. A pseudoparticle can be 
any localized gauge field configuration. For exam- 
ple we also employ pseudoparticles with a limited 
range of interaction or without topological charge. 

• In addition to a color orientation matrix we assign 
to each pseudoparticle a variable amplitude. In this 
way pseudoparticles are also able to model small 
quantum fluctuations. 

Note that the pseudoparticle approach is a general 
technique, which is in no way restricted to SU(2) Yang- 
Mills theory With minor modifications it can also be 
applied to other quantum field theories, e.g. quantum 
mechanics (c.f. 0). 



where rfy v = e bflv + 6 bfl S QL , - d b „6 Qfl and 

rftv = e b/iv - Sb^Sov + SbvSop- Each pseudoparticle has 
an index i, an amplitude A(i) G M, a color orientation 
matrix C (i) £ SO(3), a position z{i) 6 R 4 and a size 
A 6 R + . When considering a single pseudoparticle, a 
color orientation matrix is equivalent to a global gauge 
transformation. Since such a global gauge transformation 
can be specified by an element of SU(2), for which S 3 is 
a suitable parameter space, it can be expressed in terms 
of (co(i), . . • , c 3 (>')) E S 3 , i.e. c (i) 2 + c(i) 2 = 1: 

C ab (i) = 6 ab (co(i) 2 -c(i) 2 ) + 2c a (i)c b (i) + 

2e abc c (t)c c (i) (10) 

(c.f. Appendix |A"]). Note that Euclidean Lorentz transfor- 
mations are equivalent to color rotations, when consider- 
ing a single instanton or antiinstanton [1 01 ] , while a single 
akyron is form invariant under such Lorentz transforma- 
tions. Therefore Euclidean Lorentz transformations have 
been included by considering arbitrary color orientations. 

Setting A(i) = 2 in {7J yields an instanton in regular 
gauge Although for A(i) ^ 2 such pseudoparticles 
are not actually instantons, we nevertheless refer to them 
by that term. For A(i) = 2 the action of an instanton 
is S = 8ir 2 /g 2 , otherwise it is S = oo. The topological 
charge is given by 

Q = \A(i) 2 (3-A(i)). (11) 

With exception of a sign reversal in pip the same applies 
for antiinstantons (0. 

A single akyron [351 ] is a pure gauge, i.e. S = and 
Q = 0. Note that for linear superpositions of akyrons 
S ^ in general, due to the non-Abelian nature of SU(2). 
However, any such superposition has vanishing topologi- 
cal charge density (c.f. Appendix [B]) . 

A common and essential property of instantons, antiin- 
stantons and akyrons is their long range nature. For large 
| a; | the corresponding gauge fields decrease like l/\x\. As 
a consequence these pseudoparticles have the ability to 
interact over large distances. 



C. The standard choice of pseudoparticles: 
"instantons", "antiinstantons" and akyrons 

For the major part of this work we consider the follow- 
ing pseudoparticles: 

A a Jx) = A{z)C ab (i)a b wton {x-z{i)) , 



p, instanton 
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x 2 + A 2 ' 



(7) 



(8) 



(9) 



Why this particular choice of pseudoparticles? 

An important reason for considering pseudoparticles 
([7]) and (JSJ) is their similarity to regular gauge instantons 
and merons, which are known to exhibit confinement [3, 

We additionally include akyrons © so that the gauge 
field has both a transverse part and a longitudinal part 
(superpositions of instantons and antiinstantons form 
transverse gauge fields, whereas superpositions of aky- 
rons form longitudinal gauge fields; c.f. Appendix |C|) . 
Furthermore, one can show that in the continuum limit, 
i.e. in the limit of infinitely many pseudoparticles, instan- 
tons (or antiinstantons) and akyrons form a basis of all 
gauge field configurations (c.f. Appendix [Dj) . 
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Finally, numerical studies with different types of pscu- 
doparticles have shown that observables are not very 
sensitive to moderate changes in the definition of the 
pscudoparticles. For example replacing l/(a; 2 + A 2 ) by 
(1/A 2 ) exp(-a; 2 /2A 2 ) in © to © hardly affects numeri- 
cal results, if A is sufficiently large (c.f. section llV A2j) . ft 
seems that results in the pseudoparticlc approach mainly 
depend on certain "global pseudoparticlc properties" , 
like their ability to interact over sufficiently large dis- 
tances or whether they carry topological charge or not. 

D. Pseudoparticle ensembles 

We put N pseudoparticles with randomly and uni- 
formly chosen positions inside a hyperspherical spacetime 
volume (c.f. Fig. [l}. In the following we denote the ra- 
dius of this spacetime hypersphere by r spaC otimo and its 
volume by K, pace time- These quantities are related accord- 
ing to 14p aco ti m c = (7r 2 / 2 ) r s P acetime- The pseudoparticle 
density is given by n — -/V/14pacetimc- For most calcu- 
lations we consider N = 400 pscudoparticles © to © 
with size A = 0.5 and pseudoparticle density n = 1.0. 
This amounts to r space time = 3.0. 

The gauge field is a linear superposition of instantons, 
antiinstantons and akyrons: 

A-(x) = ^^)C Qb ( l )a^ instanton (x-z(*)) + 

i 

(j) a ^,antiinstanton( X — 

3 

J2A(k)C ab (k)al akymn (x~z(k)) (12) 

k 

(the indices i, j and k assume values from different 
ranges). In accordance with the ratio of transverse and 




boundary of spacetime 



FIG. 1: a 4-dimensional spacetime hypersphere of radius 
rspacetime with N = 20 pseudoparticles (pseudoparticles are 
symbolized by dark gray circles). 



longitudinal gauge field components, which is 3 : 1, we 

choOSe ^Vjnstanton ■ -^antiinstanton • -^akyron 3.3.2 

(Mnstanton, ^antiinstanton and iV akyron arc the correspond- 
ing pseudoparticlc numbers). Although in the limit of 
infinitely many pseudoparticles instantons by themselves 
form a basis of all transverse gauge fields, i.e. applying in- 
stantons and antiinstantons would yield an overcomplete 
basis, we consider an equal number of instantons and 
antiinstantons, when using a finite number of pscudopar- 
ticles. Our ensembles are then symmetric with respect 
to topological charge (c.f. (fTTjl h 

We define the ensemble average of a quantity O by 

(o) = || (^JdA(i)dC(i)^0[A(A(i),C(i))] 

e -S[A(A(i),C(i))] 

Z = J (j[dA(i)dC{i^je~ s ^ A ^ c ^ (13) 
S[A{A{i),C(i))] = f d 4 xs{A(A(i),C(i))) (14) 

" Vgpacctimc 

(dC(i) is the invariant integration measure on iS 3 ), i.e. 
the integration over all gauge field configurations in the 
path integral is replaced by an integration over the ampli- 
tudes and color orientations of the pscudoparticles. The 
quantity O, the action S and the action density s can be 
expressed in terms of A(i) and C(i) via (p~2|) . Only the 
action inside the spacetime hypersphere is considered for 
such a "path integral" . 

To eliminate the dependence of the ensemble aver- 
age (O) on the randomly chosen pseudoparticle positions 
z(i), we independently calculate (O) for many different 
sets of positions and take the average. 

Instead of integrating over amplitudes and color orien- 
tations as in p3p one could also integrate over positions 
or only over amplitudes. However, we have found that 
numerical results of such computations are very similar 
(c.f. |,ii). 

Note that our pseudoparticle ensembles are Lorentz 
invariant with respect to the center of the spacetime hy- 
persphere as well as globally gauge invariant. 

Numerical realization of pseudoparticle ensembles 

Ensemble averages (O) in the pseudoparticle approach 
are given by multidimensional integrals (c.f. (|13[1 ). which 
can be computed by Monte-Carlo methods. We have ap- 
plied the Metropolis algorithm (c.f. e.g. [I2|)- In a single 
Metropolis step the pseudoparticles are updated one by 
one in fixed order. When updating a pseudoparticle, its 
amplitude and color orientation are randomly changed. 
These changes are always accepted if they reduce the ac- 
tion. Otherwise they may be rejected depending on the 
outcome of a stochastic experiment. The action inside 
the spacetime hypersphere (|14p is approximated by stan- 
dard Monte-Carlo sampling. 
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To get rid of unwanted boundary effects, samples of 
physically meaningful quantities are always taken inside 
a hypersphcrical spacctimc region of radius 
^boundary < Spacctimc (c.f. Fig. [J). We determine 
^boundary by considering plots of the average action den- 
sity (s) against the distance to the center of the spacetime 
hypersphere: ^boundary i s assigned that distance, where 
(s) starts to deviate significantly from its constant value 
near the center of the spacetime hypersphere. In the ex- 
ample shown in Fig. [5] we have chosen ^boundary = 1-8. 

Numerical calculations have shown that pseudoparti- 
cle results are pretty stable, when applying between 100 
and 800 pseudoparticles (c.f. Q). Using significantly less 
than 100 pseudoparticles is usually not possible, because 
to extract physically meaningful results one requires a 
sufficiently large spacetime region, where boundary ef- 
fects are negligible. 



<s>(distance) - N = 400, n = 1 .0, X = 0.5, g = 4.0 
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FIG. 2: N = 400, n = 1.0, A = 0.5, g = 4.0. (s) plotted 
against the distance to the center of the spacetime hyper- 
sphere. 



where O and a are the numerical dimcnsionlcss values 
of O and a at coupling constant g and cr p h ys icai is the 
value of the string tension in physical units (through- 
out this work a p h ys icai = 4.2/fm 2 ). A crucial property 
of any trustworthy numerical method is that dimension- 
less ratios of dimensionful quantities, e.g. x 1 ^ 4 / " 1 ^ 2 or 
Icriticai/c 1 / 2 , do not depend on the coupling constant, 
i.e. that dimensionful quantities scale consistently. This 
has been observed both in lattice calculations and in the 
pseudoparticle approach (c.f. section [Hi) . 

In lattice gauge theory the ultraviolet regulator is the 
lattice spacing a. Since different values for g yield dif- 
ferent values for tr, the lattice spacing in physical units 
^physical can be adjusted by choosing appropriate values 
for the coupling constant (replace ©physical by a p h ys icai 
and set L = 1 and a = 1 in ([15)0 . In the pseudopar- 
ticle approach the minimum size of ultraviolet fluctua- 
tions is determined by the average pseudoparticle dis- 
tance d = I/71 1 / 4 and the pseudoparticle size A, i.e. there 
are two ultraviolet regulators. A variation of the cou- 
pling constant g has a similar effect on these regulators 
as it has in lattice calculations on the lattice spacing, 
that is ^physical and Aphysicai are increasing functions of g 
(c.f. Fig. [3J the scale has been set by the string tension, 
which we have obtained via generalized Creutz ratios (c.f. 
section Iffl A lj) h 

In contrast to lattice calculations we additionally have 
to specify the ratio X/d. For our standard choice of pa- 
rameters, d = 1.0 and A = 0.5, the ultraviolet regulators 
are of the same order of magnitude. Numerical calcu- 
lations have shown that physical quantities are pretty 
stable with respect to a variation of X/d in the range 
0.2... 1.1 and that there is qualitative agreement with 
lattice results for suitably chosen pseudoparticles (c.f. 
section IW A lj) . 

When considering the limit of infinitely many pseu- 
doparticles, any gauge field configuration can be repre- 
sented by a suitable superposition of instantons and aky- 



E. The pseudoparticle approach compared to 
lattice gauge theory 

In many aspects the pseudoparticle approach is simi- 
lar to lattice gauge theory. Both approaches approximate 
path integrals by integrating over a finite number of de- 
grees of freedom. In lattice gauge theory the integration 
is performed over a fixed number of link variables, while 
in the pseudoparticle approach there is a fixed number 
of pseudoparticles with amplitudes and color orientations 
as degrees of freedom. 

In both cases the scale can be set by identifying the nu- 
merical value of any dimensionful quantity with its phys- 
ical/experimental value. When using the string tension 
<7, any quantity O with dimension (length)^ is given in 
physical units by 



O 



physical 



^"physical 



L/2 



o. 



(15) 



0.3 

E 

= 0.25 
f 0.2 
1 0.15 

f 0.1 

a. 

0.05 




(g) and ^„ hvsica i(g) in fm - N = 400, n = 1.0, X = 0.5 



.••**" 



.....r 



i 



average pseudoparticle distance d physica | in fm 
pseudoparticle size ^-physical in fm 

1 2 3 4 5 

coupling constant g 



FIG. 3: N = 400, n = 1.0, A = 0.5. d phys icai in fm and 
Aphysicai in fm plotted against g. 
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rons (c.f. section [D]), i.e. in this limit our pseudopartiele 
ensembles are gauge invariant. However, unlike in lattice 
gauge theory we only have approximate gauge invariancc 
when working with a finite number of degrees of freedom. 
Whether this approximate gauge invariance is a problem 
can be tested by numerical calculations. A suitable way 
of doing that is to fix the gauge and to perform com- 
putations similar to those without gauge fixing. In the 
pseudopartiele approach a particularly convenient gauge 
is Lorentz gauge, i.e. d^A^ = 0, which amounts to con- 
sidering only instantons and antiinstantons but no aky- 
rons. The qualitative agreement between these results 
and results without gauge fixing (c.f. section HVBp in- 
dicates that in the pseudopartiele approach the lack of 
exact gauge invariance is not a problem. 

In contrast to lattice gauge theory we always work in 
the continuum, i.e. we do not discretize spacetime. This 
might be an advantage when calculating glueball masses 
from correlation functions, because due to exact rota- 
tional symmetry it is possible to project out states with a 
well defined angular momentum. Furthermore, the pseu- 
dopartiele approach has a high potential when dealing 
with fermionic problems, because there is no fermion dou- 
bling in such a continuum formulation. 

Another difference lies in the number of degrees of free- 
dom, which is significantly smaller in the pseudopartiele 
approach than it is in typical lattice calculations (1,600 
degrees of freedom, when using 400 pseudoparticles, com- 
pared to e.g. 786, 432 degrees of freedom on a "small" 
16 4 -lattice). Therefore, when pseudopartiele results are 
in agreement with results from lattice calculations, the 
pseudoparticles are chosen such that they represent es- 
sential degrees of freedom of SU(2) Yang-Mills theory. To 
put it another way, by applying different types of pseu- 
doparticles one can determine the effect of certain classes 
of gauge field configurations on certain observables, and 
by doing that try to find out essential properties of con- 
fining gauge field configurations (c.f. section ITV]). Lattice 
calculations, on the other hand, are much faster, i.e. they 
can cope with a significantly larger number of degrees of 
freedom in comparable computation time. Therefore, lat- 
tice calculations are certainly better suited for producing 
high quality numerical results. 



III. CALCULATING OBSERVABLES 

A. The static quark antiquark potential at zero 
temperature 

The common tool for studying the potential of a pair 
of infinitely heavy quarks at zero temperature are Wilson 
loops. A Wilson loop W z is defined by 



W Z [A) = ^Tr( 



(16) 



where z is a closed spacetime curve and P denotes path 
ordering. Rectangular Wilson loops with spatial exten- 



sion R and temporal extension T are denoted by W^t) ■ 
It is well known that the potential of a static quark 
antiquark pair Vq q with separation R can be related to 
ensemble averages of rectangular Wilson loops according 
to 



V m {R) 



lim — In 



ir. 



(R,T) 



(17) 



(c.f. e.g. Li). _ 

In the following we assume that for large separations 
the static quark antiquark potential can be parameter- 
ized by 



V m (R) 



V 



R 



(jR. 



(18) 



This ansatz is based on the bosonic string picture [HI, 03] 
and on various numerical results from lattice calculations 
(c.f. e.g. [ID, [lj| Qj})- There are three parameters: 

• Vq is a constant shift of the potential without phys- 
ical relevance. 

• a is the coefficient in front of the attractive \/R- 
correction of the potential for large separations. 
The bosonic string pictu re p redicts 

^string = 7r/12 « 0.26 [U [3| . In lattice calcu- 
lations Qiattice = 0.22 ... 0.32 has been observed 
[IIIHzl- We refer to a as Coulomb coefficient. 

• The string tension a characterizes the force be- 
tween a static quark and a static antiquark at large 
separations. Lattice calculations yield a positive 
value of the string tension, which implies the pres- 
ence of confinement (c.f. e.g. 0,0, 13 [ill). Fur- 
thermore, a is a monotonically increasing function 
of the coupling constant g, i.e. when the scale is 
set by the string tension, the extension of the lat- 
tice in physical units can be adjusted by choosing 
appropriate values for g. 



1. Calculating a and a via generalized Creutz ratios 

We determine the string tension a and the Coulomb 
coefficient a by g eneralizing the well known method of 
Creutz ratios [19l ]. 

Generalized Creutz ratios are based on guessing the 
functional dependence of ensemble averages of rectangu- 
lar Wilson loops. A possible candidate is given by 



- In ( W, 



(R,T) 



= V [R + T - a 



R T\ 
T + R^ +fj 



crRT 



(19) 



[l6j ]. It is a simple and plausible choice, which is consis- 
tent with numerical results from lattice calculations. It 
also fulfills the following necessary requirements: 

• {W(r.t)) = {W(t,r)) (in Euclidean spacetime there 
is no difference between space and time). 
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• lim T ^oo(-ln(W (il) T)}) = V m (R)T (c.f. $T7§ and 
(USD- 
Note that (HU) is only valid for R,T>d 7 A. Wilson loops, 
which arc significantly smaller than the two ultraviolet 
regulators d and A, are subject to strong cutoff effects. 

When applying generalized Creutz ratios, the starting 
point is a set of ensemble averages of rectangular Wilson 
loops, 



W(Hi,Ti)) ' ( W (R2,t 2 ; 



(Rn,T„ 



}, (20) 



with at least a couple of different ratios Ri/Ti. 
A generalized Creutz ratio Tx is defined by 



Tx (Rii > Ti t , . . . , Ri 4 , Tj 4 ) 

Cl,X 



= (w, 



tv, 



(Ri 3 ,Ti 3 ) 



(Hi 4 ,T< 4 ) 



(21) 



with weights c Xt x = ci,x (Rh , , . . . , R ii , T i4 ), . . . , 
C4,a" = ca,x{Rii i Ti x , . . . , Ri 4 , Tj 4 ), which will be specified 
below (X denotes any of the constants Vq, a, (3 or a). 
Inserting the Wilson loop ansatz (p~9|) yields 



In ( I\y (-Rii , Ti ± , . . . , i?i 4 , T i4 ) 
V (c ltX (R il +T h 



C4,X ( -Ri 4 + T i4 



C4,X 



Rii Ti\ 

■Mi -Rii 
-Ri 4 ^i 4 
-M 4 -Ri 4 



''^<- .v • ... • '•; A J + 

o ( ci,x (RiiTi^J + ■ • ■ + c 4i x ( RuT u 



(22) 



By solving a linear system the weights ci^x, ■ ■ ■ , c^,x are 
chosen such that (l22l) reduces to 



(23) 



-In rx^,^,...,^,^) = X 



To determine Vq, a, (3 and c, we consider all four ele- 
ment subsets of (|20|) . For every subset we calculate Vb, 
a, /? and a via (|2"3"|) . Although there will be certain fluc- 
tuations due to systematical and statistical errors, most 
of these estimates should be pretty similar. 

We have found that estimates for Vb, a, (3 and a with 
large x + . . . + c\ x exhibit significantly stronger fluc- 
tuations than estimates with small cf x + . . . + c\ x (c.f. 
Fig. [4J. Therefore, we sort the estimates according to 
cf x + . . . + c\ x and only keep the "smaller half" . The 
average of these estimates is the final result for Vb, a, (3 
or a . 



Results: N = 400, n = 1.0, A = 0.5 

We have computed ensemble averages of rectangular 
Wilson loops (W ( r.t)) for R,T & {7a , 8a , . . . , 12a}, 
a = 0.21, ?*boundary = 1-8, and have determined all pos- 
sible generalized Creutz ratios. There are around 6, 000 
generalized Creutz ratios for each of the quantities Vq, a, 
(3 and a. 

In Fig. 2] the corresponding estimates for the string 
tension a and the Coulomb coefficient a are plotted for 
g — 5.0, sorted according to c\ x + . . . + c\ x in increasing 
order from left to right (for the sake of clarity only every 
tenth estimate is shown). Both plots demonstrate that 
estimates on the left exhibit a smaller variance than es- 
timates on the right. This indicates that estimates with 



small x 
with large c\ x 



"4,X 



are more reliable than estimates 



-a,x- 



To extract numerical val- 



ues for a and a, we have calculated the average and the 
standard deviation from the "smaller half" . The results 
are er = 0.342 ± 0.016 and a = 0.232 ± 0.022. 

By proceeding in the same way we have obtained nu- 
merical values for a and a for different coupling constants 
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FIG. 4: N = 400, n = 1.0, A = 0.5, g = 5.0. a) Estimates for 
a sorted according to c\ x + . . . + c| x . b) Estimates for a 
sorted according to cf x + • ■ • + c l x ■ 
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g€ {2.0, 2.5, ... , 5.5}. 

In Fig. [5^ the string tension a is plotted against the 
coupling constant g. It is positive and an increasing func- 
tion of g. When the scale is set by the string tension, the 
size of the spacetime hypersphere in physical units can be 
adjusted by choosing appropriate values for the coupling 



We compute ensemble averages of rectangular Wilson 
loops (W^t)) for fixed R but different T to obtain 
a curve — ln(M / (jj =consta nt,T)}- Examples are shown in 
Fig. [B^i. According to such a curve will exhibit a 
linear behavior for sufficiently large T. From the slope, 
which we obtain by fitting a straight line, we can read 



constant. For a 



physical 



4.2/fm its diameter ranges off Vc,q(-R). Iterating this procedure for a number of dif- 



from approximately 0.81 fm at g = 2.0 to 1.92 fm at 
g = 5.5. 

Fig. E}j shows the Coulomb coefficient a as a function 
of the coupling constant g. It increases from a = 0.04 at 
g = 2.0 to a = 0.25 at g = 5.5. For large g the value 
of a is in qualitative agreement with the prediction from 
the bosonic string picture, String =tt/12« 0.26 [TEEl, 
and with results from lattice calculations, 
lattice = 0.22... 0.32 [IE 03. 

a) o(g) - N = 400, n = 1.0, X = 0.5 




12 3 4 

coupling constant g 



N = 400, n = 1.0, X = 0.5 



e 

o 

o 
o 




2 3 4 

coupling constant g 



FIG. 5: N = 400, n = 1.0, A = 0.5. a) a plotted against 
b) a plotted against g. 



2. Calculating the static quark antiquark potential 

The starting point to calculate the static quark anti- 
quark potential is (|17| . which we write as 



In < W, 



(fl,T) 



(24) 



ferent R yields an approximation of the quark antiquark 
potential. 



Results: N = 400, n = 1.0, A = 0.5, g = 5.0 

We have computed ensemble averages of rectangular 
Wilson loops (W(h,t)> for R,T 6 {a, 2a , ... , 12a}, 

a = 0.21, rboundary = 1-8. 

Fig. [Si shows - ln(W (J?=constant!T) ) as a function of T 
for different R 6 {a , 2a , . . . , 12a}. To determine the 
slope of these curves for large T, we have fitted straight 
lines to the data points at T € {9a , 10a , 11a , 12a}, as 
indicated by the dashed straight lines. 

The corresponding potential as a function of the quark 



a) -ln<W (aT) >| R = constant (T) - N = 400,n = 1.0,A, = 0.5,g = 5.0 
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b) V (R) - N = 400, n = 1.0, X = 0.5, g = 5.0 



0.8 - 



0.6 



least squares fit with V - a/R + aR 
cutoff effects 



0.5 



1 1.5 2 

quark antiquark separation R 



2.5 



FIG. 6: JV = 400, n = 1.0, A = 0.5, g = 5.0. 

a) — ln(W / ( H=constant T )) for different R plotted against T 

(a — 0.21). b) Kjq plotted against the separation of the 

quarks. 
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antiquark separation is plotted in Fig. \6jp. For large sepa- 
rations it clearly exhibits a linear behavior. To obtain nu- 
merical values for the string tension a and the Coulomb 
coefficient a, we have performed a least squares fit of 
the potential parameterization (|18p to the data points 
shown in Fig. [BJd. Only data points with R > 1.3 > d, A 
have been considered, because cutoff effects arc expected 
to render the potential unphysical for small separations. 
Within statistical errors the results, a = 0.376±0.026 and 
a = 0.283 ± 0.093, are in agreement with the results ob- 
tained by generalized Creutz ratios (c.f. section [Til A f[) . 
The fit is also shown in Fig. [Sb- 

Note that in contrast to generalized Creutz ratios we 
have not made any assumption about the functional de- 
pendence of ensemble averages of Wilson loops to calcu- 
late the quark antiquark potential. Therefore, the agree- 
ment of the results for a and a with results obtained by 
generalized Creutz ratios shows again the consistency of 
the Wilson loop ansatz (fT9|) and Monte-Carlo data for 

(W(fl,T))- 



B. The topological susceptibility 

In order to produce quantitative results involving the 
string tension, we need other dimensionful quantities so 
that we can consider dimcnsionlcss ratios. One such 
quantity, which has been studied extensively on the lat- 
tice (c.f. e.g. [13, HH, HI, HH), is the topological suscepti- 
bility x- The topological susceptibility is closely related 
to the mass of the rf meson [2J|. It is defined by 



X 



(25) 



where Qv is the topological charge inside the spacetime 
volume V. 

In our numerical calculations we approximate the limit 
V — > oo by a finite volume, the hyperspherical spacetime 
region with radius ^boundary (c.f. section III D I) . A number 
of computations has shown that x is stable with respect 
to a a variation of that spacetime volume. 



Results: N = 400, n = 1.0, A = 0.5 



Fig. [7| shows the dimensionless ratio x^ 4 / ^ 2 as a 
function of the coupling constant g (^boundary = 1-8, 
<j has been obtained by generalized Creutz ratios (c.f. 
Fig-Ek))- As expected, x 1 ^ 4 / " 1 ^ 2 is nearly independent 
of g, i.e. the string tension and the topological suscep- 
tibility exhibit consistent scaling behaviors with respect 
to the coupling constant. This success strongly indicates 
that the pseudoparticle approach has the potential to re- 
produce correct Yang-Mills physics. The range of values, 
X 1 / 4 /c 1 / 2 = 0.33 . . .0.35, is of the right order of magni- 
tude compared to the lattice result 



0.5 ■ 



0.3 



0.2 



0.1 - 



X 1/4 to 1/2 (g) - N = 400, n = 1.0, X = 0.5 
lattice result: jc 1/4 /a 1 ' 2 = 0.486 ± 0.010 




1 2 3 4 5 6 

coupling constant g 

FIG. 7: N = 400, n = 1.0, A = 0.5. x 1/4 /^ 1/2 plotted against 
9- 



C. The critical temperature of the confinement 
deconflnement phase transition and the finite 
temperature quark antiquark potential 

1. The pseudoparticle approach in finite temperature SU(2) 
Yang-Mills theory 

For finite temperature calculations we have to con- 
sider periodic gauge field configurations in xo-direction. 
Therefore we need a method to generate periodic building 
blocks. In the case of singular gauge instantons this has 
already been done in form of calorons [I HI, which are 
periodic solutions to the classical Yang-Mills equations 
of motion. However, there is no straightforward gener- 
alization to arbitrary pseudoparticlcs. Therefore, we in- 
troduce a different method, which resorts to a blending 
technique from computer aided geometric design 27 1. 

At first, we define blending functions B and B, which 
form smooth, C 1 -continuous connections between at 
A = and 1 at A = 1 and vice versa: 



B(X) = -2A 3 + 3A 2 



B(\) 



l-B(X). (26) 



To make a non-periodic pseudoparticle a° with its center 
at the origin periodic in xo-direction, we multiply "both 
ends" with the blending functions and add the results 
(c.f. Fig.©: 



7° (r) = 

/^periodic V / 

a£0co,x) if 



B{X)a a ^ {x -j9,x) + £(A)< (*o , x) 



A = 



xo-ifl- b)/2 



To evaluate 



/^.periodic 



at Xq 



(27) 

-{/3-b)/2,(l3 + b)/2} one 



just has to combine (|2~7| and 



( X V4 /(J l/ 2) 



lattii 



0.486 ±0.010 \2t 



a 



u, periodic 



(x ,x) = a a iodic (x Q + n/3, x) 



(28) 
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blending region b blending region b 



FIG. 8: the basic principle of the blending method. 



with a suitably chosen integer n. b is the width of the 
blending region. Throughout this work we have chosen 
b = 0.3 x (3. 



The spacetime region 

At zero temperature we consider a spacetime region, 
which is the interior of a 4-dimensional hypersphere (c.f. 
Fig. [1]). At finite temperature that spacetime region is 
replaced by a spacetime with a periodic time direction 
of extension (3 and a spatial part, which is the interior 
of an ordinary 3-dimcnsional sphere of radius r spa co (c.f. 
Fig. El). 

As in the zero temperature case we have to assure that 
samples of physically meaningful quantities are always 
taken inside a spacetime region, where boundary effects 
can be neglected. The spatial part of such a region is 
the interior of a sphere of radius ?*b oun dary! whereas the 
time direction is not restricted (the light gray region in 
Fig.©. 



regions of strong region of negligible 




boundary of space 
(3-dimensional sphere) 

FIG. 9: the spacetime region at hnite temperature. 



The pseudoparticle density 

At hnite temperature we use the same pseudoparticle 
density as for the corresponding zero temperature cal- 
culation (a zero temperature calculation with identical 
parameters is necessary to obtain a numerical value for 
the string tension, which is used to set the scale). This 
is a reasonable choice, because in the limit (3 — » oo, 
^spaco — > oo and r spaco ti mc — > oo hnite temperature en- 
sembles and zero temperature ensembles become identi- 
cal. Furthermore, there is a close analogy to lattice calcu- 
lations. To set the scale, a zero temperature calculation 
with the same number of lattice sites in all four spacetime 
directions is carried out. For hnite temperature calcula- 
tions the number of lattice sites in the temporal direction 
is significantly reduced. However, the density of lattice 
sites is still the same, i.e. there is a smaller number of 
link variables in a smaller spacetime volume. 

The pseudoparticle approach offers two possibilities to 
change the temperature: 

• Varying (3 while g is kept constant changes the value 
of the temperature directly. 

• Varying g while (3 is kept constant changes the nu- 
merical value of the string tension a. This leads to 
a different extension of the periodic time direction 
in physical units and, therefore, changes the value 
of the temperature. 

Both methods yield consistent results (c.f. Fig. [TUbh 



2. The critical temperature of the confinement 
deconfinement phase transition 

Whereas at low temperatures static quarks cannot be 
separated from each other, at high temperatures a de- 
confining phase is expected. The temperature, at which 
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the corresponding phase transition takes place, is called 
critical temperature and denoted by T C riticai- 

^critical can be determined from ensemble averages of 
Polyakov loops. A Polyakov loop is a Wilson loop around 
the periodic time xq: 



P< exp [i <b dz A (z) 



(29) 



Due to spatial translational invariance {L z )p is z-indc- 
pendent. Therefore, from a numerical point of view it is 
convenient to consider spatial averages of Polyakov loops. 
We define 



L 



d 3 zL z [A\ 



(30) 



(L)p can be considered as an order parameter indicating 
whether there is confinement or not (c.f. e.g. [l2l. l28j|): 



= 



L) / 





confinement 



deconfinement . 



(31) 
(32) 



This criterion is closely related to center symmetry, which 
is spontaneously broken in the dcconfincment phase. 

In the center symmetric phase of SU(2) Yang-Mills the- 
ory for every field configuration A a ^ + there is a gauge 
equivalent field configuration A a ^~ with 
L Z [A + ] = -L Z [A~] (A^ + and A^~ are related by center 
symmetry, i.e. they are connected by a singular gauge 
transformation) . Since both field configurations have the 
same action, i.e. 5[^4 + ] = (L)p = follows imme- 

diately. 

The spontaneous breakdown of center symmetry at 
T = T critica i comes along with a splitting of the Hilbcrt 
space of states in two independent spaces: 
TL — ► H + © H~ . The same applies for the set of field 
configurations considered in the path integral: 
A — > A + © A~ . The ensemble average of the Polyakov 
loop depends on which Hilbert space was chosen dur- 
ing the spontaneous breakdown of center symmetry: 
( l )/3,h+ = +~ l and (L) PM - = -I. 

In the broken phase two field configurations, which are 
related by center symmetry, cannot be connected con- 
tinuously by a set of field configurations of finite action. 
This implies that during a Monte-Carlo simulation only 
field configurations either from A + or from A~ are gener- 
ated, assuming an infinite system and a local and contin- 
uous update mechanism. In numerical calculations these 
assumptions are only approximately fulfilled. Neverthe- 
less, in lattice Monte-Carlo simulations (c.f. e.g. [H, [29|) 
it has been observed that there are long sequences of 
steps, where only field configurations corresponding to 
one of the two Hilbert spaces are generated. 

Since our pseudoparticle ensembles are only approx- 
imately center symmetric, there is a smooth transition 
between the two phases. (L)p m can be observed well 
below the critical temperature. For (3 « /^critical the en- 
semble average of the Polyakov loop quickly rises to a 



non-zero value. For high temperatures {L)p ~ 1 (c.f. 
Fig. 1 10b). Therefore, we define the critical temperature 
^critical (or equivalently its inverse /^critical) to be that 
temperature, where the ensemble average of the Polyakov 
loop crosses a certain value £, i.e. 



L 



(33) 



{L)p in the pseudoparticle approach versus {L)p in lattice 
calculations 



The main difference between our results and lattice re- 
sults is that even in the deconfinement phase we never 
have observed (L)p < 0. We conclude that in the pseu- 
doparticle approach the low action field configurations 
corresponding to Hr are underrepresented, i.e. there are 
more low action field configurations corresponding to H + 
than to This bias always forces the system to chose 
the Hilbcrt space T~l+, when center symmetry is sponta- 
neously broken. 

The bias can be explained by the following qualitative 
argument: for field configurations close to zero L z w 1 
(c.f. (f2"9"| ). Therefore, we expect these field configurations 
to be elements of A + . Furthermore, all of these field con- 
figurations are low action field configurations, which con- 
tribute significantly to the path integral. The set of these 
field configurations is denoted by Aa^o- Loosely speak- 
ing, all other low action field configurations are large 
enough so that L z can pick up exponents, which are sig- 
nificantly different from zero. They are cither elements 
of A + or of A~ . The set of these field configurations 
is denoted by Aa^o- However, there is numerical evi- 
dence that in the pseudoparticle approach with around 
400 pscudoparticlcs most low action field configurations 
have small gauge fields, i.e. are elements of Aa^o- This 
is closely related to the fact that our ensembles are only 
approximately gauge invariant: a gauge transformed field 
configuration from Aaxio, which is a field configuration 
in Aa^>o, can only be approximated, when using around 
400 pseudoparticles; such an approximation usually has 
a higher action. Therefore, there are proportionally more 
field configurations in Aa^o than in Aa^>o- This is man- 
ifested in form of a bias. 



Results: 



= 3.00, n = 1.0, A = 0.5 



In Fig. [TUk the ensemble average of the Polyakov loop 
(L}/3 is plotted against the temperature T for different 
coupling constants (r spacc = 3.00, r bo undary = 1.8). The 
dashed straight lines correspond to a determination of 
^critical with £ = 0.5 (other values for £ yield qualitatively 
identical results). 

Fig. [10b shows that the dimensionlcss ratio 
Tcriticai/c 1 / 2 is nearly independent of the coupling 
constant g (a has been obtained by generalized Creutz 
ratios (c.f. Fig. E^)), i.e. the string tension and the 
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a) (L)p(T) - r space = 3.0,n = 1.0, X = 0.5 
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coupling constant g 



methods to adjust the temperature, varying f3 or varying 
g. Fig. [TUb shows (L)^ as a function of T/cr 1 / 2 for differ- 
ent coupling constants g. The fact that all sample points 
scale to a single curve demonstrates that both methods 
yield consistent results. 

Fig. [TT] shows the evolution of during a single 
Monte-Carlo simulation (1,000 Monte-Carlo steps) for 
g = 5.0 and different T/T crit icai (^critical = 0.38). These 
plots demonstrate that in the confinement phase 
(T < Td-iticai) L z assumes approximately an equal num- 
ber of positive and negative values, whereas in the decon- 
finement phase (T > T cr iti ca i) the values of L z are mainly 
positive. 




c) <L> p (T/a 1 ' 2 ) - r space = 3.0, n = 1 .0, X = 0.5 




T/T c 



critical 



= 1.50 



FIG. 11: r spacc = 3.0. 
different ratios T/T cr iticai 
Monte-Carlo step. 



n = 1.0, A = 0.5, f 

plotted against the number of the 



5.0. L 3 for 



3. The static quark antiquark potential at finite temperature 

It is well known that the static quark antiquark poten- 
tial at finite temperature can be obtained from Polyakov 
loop correlation functions: 



1 



FIG. 10: r spaco = 3.0, n = 1.0, A = 0.5. a) {L) p for different g V m( R > I 3 ) ~ p ln \ L o L 



plotted against T. b) Tcriticai/c plotted against g. c) (L)g 
for different g plotted against T/a 1 ^ 2 . 



critical temperature scale consistently with respect to 
g. We would like to stress again that such a scaling 
behavior, although mandatory for any trustworthy 
numerical method, is far from obvious. The range 
of values, Tcriticai/c 1 / 2 = 0.54 . . . 0.65, is of the right 
order of magnitude compared to the lattice result 

(^critical/ O ' ) lattice — 

0.694 ±0.018 \M- 
As we have explained in section IIII C 11 there are two 







= R 



(34) 



(c.f. e.g. 0,121). 

However, center symmetry is only approximately re- 
alized in the pseudoparticle approach, i.e. even in the 
confinement phase {L)p ^ 0. Consequently, any finite 
temperature potential calculated according to (|34[) ap- 
proaches a constant value for large quark antiquark sep- 
arations: 

2 



lim V m (R,p) 



In 



L 



(35) 



Therefore, at finite temperature one can expect to ob- 
serve confining potentials only up to intermediate sepa- 
rations. 
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Results: 



= 3.00, n = 1.0, A = 0.5, g = 5.0 



A. Pseudoparticles of different size and profile 



Fig. [T^] shows finite temperature quark antiquark po- 
tentials for T/T critical G {0.39 , 0.52 , . . . , 2.09} obtained 
from Polyakov loop correlation functions (c.f. (f3"4")) ; 
^spacc = 3.00, /"boundary = 1-27) and the corresponding 
zero temperature potential obtained from ensemble aver- 
ages of Wilson loops (c.f. section HlIA2p . At high tem- 
peratures there is clearly no confinement. However, with 
the decrease of the ensemble average of the Polyakov loop 
towards lower temperatures finite temperature potentials 
gradually approach the zero temperature potential ex- 
tracted from Wilson loops. 

Because of (L)p ^ the finite temperature potentials 
and the zero temperature potential differ at large sepa- 
rations even for temperatures significantly smaller than 
^critical- There are arguments suggesting that these po- 
tentials should rather be calculated via 



V m (R,(3) = -±ln 
\z\ = R 



(L S - (L) p ) (Ll - (L)p 



(36) 



instead of ([34]) . Indeed, this yields finite temperature 
potentials, which are identical to the zero temperature 
potential within statistical errors, a result, which is in 
agreement with results from lattice calculations (c.f. e.g. 
[30|). We plan to discuss these topics in more detail in 
an upcoming paper. 
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FIG. 12: r spacc = 3.0, n = 1.0, A = 0.5, g = 5.0. Finite tem- 
perature potentials and the corresponding zero temperature 
potential plotted against the separation of the quarks. 



IV. PROPERTIES OF CONFINING GAUGE 
FIELD CONFIGURATIONS 

The goal of this section is to identify essential proper- 
ties of gauge field configurations, which are responsible 
for confinement. To this end we compare different pseu- 
doparticle ensembles. 



1. Pseudoparticles of different size 

In the following we explore how pscudoparticlc results 
are affected by a variation of the pseudoparticle size A, 
while the pseudoparticle density n or equivalently the 
average pseudoparticle distance d = 1/n 1 / 4 is kept con- 
stant. In other words, we consider different ratios of the 
two ultraviolet regulators d and A. Note that A strongly 
affects the shape of a pseudoparticle near its center, while 
it has essentially no effect on the long range behavior of 
a pseudoparticle, which is proportional to l/\x\ (c.f. ([7]) 
to ©). 

We consider ensembles with N = 400, n = 1.0, g = 4.0 
and A G {0.20, 0.35, ... , 1.10}. 

We have obtained the string tension a via generalized 
Creutz ratios (c.f. section UlI A 1[) . Fig. [T3"h shows that 
<7 is clearly positive for all values of the pseudoparticle 
size A. Moreover, there is only a weak A-depcndence, i.e. 
a = 0.23 . . . 0.29 for A = 0.20 . . . 1.10. The implication 
is that confinement is connected to the \j\x\ long range 
behavior of the pseudoparticles, which is the same for all 



a) o(X) - N = 400, n = 1 .0, g = 4.0 
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FIG. 13: N = 400, n = 1.0, g = 4.0. a) a plotted against A. 
b ) X 1/4 /°" 1/2 plotted against A. 
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values of A. The shape of the pseudoparticlcs near their 
center, which is A-dependent, has essentially no effect on 
the string tension and, therefore, is of no relevance to 
confinement. 

The topological susceptibility \ is also fairly indepen- 
dent of the pseudoparticle size A as is the dimension- 
less ratio x 1 ^ 4 /c 1 ^ 2 (c.f. Fig. [T3"b). The range of values, 
X 1/4 /er 1/2 = 0.37 . . . 0.31 between A = 0.20 . . . 1.10, is in 
qualitative agreement with the lattice result 



(x 1/4 /- 1/2 ) 



lattice 



0.486 ±0.010 \2l 



Pseudoparticles of Gaussian localized profile and 
different size 



a) - Gaussian pseudoparticles: N = 400, n = 1 .0, g = 4.0 



c 0.3 




To learn more about the interrelation between the long 
range behavior of pseudoparticlcs and confinement, we 
study pseudoparticles of Gaussian localized profile and 
different size. Such Gaussian localized pseudoparticles 
are obtained by replacing l/(x 2 + A 2 ) appearing in the 
definitions of instantons, antiinstantons and akyrons ((J7J) 
to ©) by (l/A 2 ^-* 2 /^ 2 ). The main difference is the 
long range behavior of the resulting building blocks: in 
contrast to our standard choice of pseudoparticlcs, ([7]) 
to ©, Gaussian localized pseudoparticles have a limited 
range of interaction, which is proportional to their size 
A. 

We consider ensembles with TV = 400, n = 1.0, g = 4.0 
and A G {0.25, 0.50, ... , 1.50}. 

We have obtained numerical values for the string ten- 
sion a via generalized Creutz ratios (c.f. section Ull A 1[) . 
The results are plotted in Fig. [T4"b against the pseu- 
doparticle size A. There is no confinement for A < 0.75 
and confinement for A > 1.00. The onset of confinement 
takes place somewhere around A ~ 0.75. This is precisely 
the width at which neighboring pseudoparticles start to 
overlap significantly (this can be seen by assigning each 
pseudoparticle an appropriate volume, e.g. the volume of 
a hypersphere of radius A, (7r 2 /2)A 4 ~ 4.93 x A 4 , and 
comparing that volume with the maximum volume non- 
overlapping pseudoparticles can cover, i.e. 1/n = 1.0). 
The implication is that pseudoparticle ensembles only 
exhibit confinement if their pseudoparticles cover suffi- 
ciently large spacetime regions so that there are signifi- 
cant overlaps between neighboring pseudoparticles. This 
is reminiscent of percolation. Note that percolation phe- 
nomena have been related to confinement in a variety of 
ways, eg. percolation of center vortices [3l[ or monopole 
loops [32j in lattice gauge theory or random percolation 
of bonds or sites on three dimensional lattices [33[ . 

The A-dependcncc of the topological susceptibility \ 
is much weaker than the A-dependence of the string ten- 
sion a. Therefore, the dimensionless ratio 

x l/4 /(T l/2 ig 

dominated by a (c.f. Fig. [T4"b). The range of values for 
A > 1.00, i.e. for significantly overlapping pseudoparti- 
cles, x 1/4 /cr 1/2 = 0.41... 0.28, is of the right order of 
magnitude when compared to the lattice result 

(x l/4 /(T l/2 )iatt . cc = ^ 486 ± a()10 
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FIG. 14: Gaussian localized pseudoparticles, N = 400, n = 
1.0, g = 4.0. a) a plotted against A. b) x 1/4 /°" 1/2 plotted 
against A. c) Kjq for different pseudoparticle sizes plotted 
against the separation of the quarks. 



We have also calculated quark antiquark potentials for 
different values of A (c.f. section Ull A 2[) . The results are 
shown as functions of the quark antiquark separation in 
Fig. [T4b. For A > 1.25 the potential is clearly confining, 
whereas for A < 0.50 it is unambiguously not confining, 
a result which is in agreement with previous results from 
this section. 
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3. "Singular gauge pseudoparticles " of different size 

We have also studied ensembles of singular gauge in- 
stantons, antiinstantons and akyrons. These "singu- 
lar gauge pseudoparticles" are obtained by replacing 
l/(x 2 + \ 2 ) in © to © by X 2 /(x 2 + e 2 )(x 2 + X 2 ), where 
e is an ultraviolet regulator, which has been included due 



a) a(X) - Singular pseudoparticles: N = 400, n = 1 .0, ... 




1 1.5 
pseudoparticle size X 



to numerical reasons. 

We consider ensembles with N = 400, n = 1.0, e = 0.1, 
g = 4.0 and A E {0.25 , 0.60 , . . . , 2.00}. 

Numerical results (c.f. Fig. fl~5|) and their interpreta- 
tion are similar to those obtained with Gaussian local- 
ized pseudoparticles. Of course, this is what one would 
expect , because the gauge fields of singular gauge pseu- 
doparticles decrease much faster than those of regular 
gauge pseudoparticles, i.e. like 1/ 1 a? | 3 instead of l/\x\. 
Therefore, singular gauge pseudoparticles are quite simi- 
lar to Gaussian localized pseudoparticles with respect to 
long range interactions. 

The result that there is no confinement for ensembles 
of singular gauge pseudoparticles of small size (X/d< 0.6) 
is in agreement with what has been obtained in instan- 
ton gas and instanton liquid models, where typically 
X/d « 1/3 (c.f. e.g. [l|). Furthermore, it has been pro- 
posed that large size singular instantons can produce con- 
finement ([34| quoted in [H), which is also supported by 
our findings. 



B. Instantons and antiinstantons versus akyrons 
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FIG. 15: Singular gauge pseudoparticles, N = 400, n — 1.0, 
e = 0.1, g = 4.0. a) a plotted against A. b) x 1/4 /o- 1/2 
plotted against X. c) F q q for different pseudoparticle sizes 
plotted against the separation of the quarks. 



In the following we present evidence that confinement 
arises due to instantons and antiinstantons, whereas aky- 
rons do not produce confinement. To this end, we con- 
sider ensembles with the same number of pseudoparticles 
but different ratios iVi nst anton : Nantiinstanton : Nakyron- In 
detail, we compare the following ensembles with 
N = 400, n = 1.0, A = 0.5 and g = 4.0: 

• "Akyron ensemble": an ensemble containing 400 
akyrons (no instantons or antiinstantons). 

• "Standard ensemble" : an ensemble containing 150 
instantons, 150 antiinstantons and 100 akyrons. 

• "Instanton ensemble" : an ensemble containing 200 
instantons and 200 antiinstantons (no akyrons). 

Fig. 1161 shows the quark antiquark potential as a func- 
tion of the separation for all three ensembles (c.f. sec- 
tion [IirX2]). The standard ensemble and the instanton 
ensemble exhibit a confining potential, whereas the aky- 
ron curve indicates that there is no confinement. 

We have obtained numerical values for the string ten- 
sion a via generalized Creutz ratios (c.f. section UlI A 1|) . 
The results are Cakyron = 0.019 ±0.08, 

^standard = 0.236 ± 0.013 and (Tinstanton = 0.512 ± 0.023. 

The ratio of these values is given by 

f^akyron . fJstandard • ^"instanton ~ 1 ■ 12 . 27. 

To compare dimcnsionless ratios, we have calculated 
the topological susceptibility \ an d the critical tempera- 
ture Tcritical in the standard ensemble and in the instan- 
ton ensemble: 

( x 1 /4 /(T i/2 )standard = ,35 5 

(X 1/4 /^ 1/2 )instanton = 0.26, 
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V qq (R) - N = 400, n = 1.0.X = 0.5, g = 4.0 
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> cutoff effects 
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FIG. 16: akyron ensemble (400 akyrons), standard ensem- 
ble (150 instantons, 150 antiinstantons and 100 akyrons) and 
instanton ensemble (200 instantons and 200 antiinstantons), 
N = 400, n = 1.0, A = 0.5, g = 4.0. V q q plotted against the 
separation of the quarks. 



(Ic r itical/o' 1 ^ 2 )standard =0.61, 
(^critical/^" ^ )instanton — 0.59. 

Obviously, akyrons increase the topological susceptibil- 
ity, while they do not affect the critical temperature. 
Since (x 1//4 /cr 1 ^ 2 ) s tandard = 0.35 is closer to the lattice 
result (x 1/4 /cr 1/2 )iatticc = 0.486 ± 0.010 [H than 
(x 1 ^ 4 /fT 1 / 2 )i nst anton = 0.26, it is beneficial with respect 
to quantitative results to consider ensembles containing 
not only instantons and antiinstantons but also akyrons. 

The dimensionless ratios x 1 ^ 4 / ' 1 ^ 2 an d T cr iti C ai/o' 1 ' /2 
in the akyron ensemble arc not meaningful: the topolog- 
ical susceptibility x vanishes identically (c.f. Appendix |B|) 
and the ensemble average of the Polyakov loop (L)p is 
close to 1 even for very large values of (3, i.e. there is 
no sign of a confinement deconfinement phase transition. 
From that and the fact that the string tension is more 
than ten times smaller in the akyron ensemble than in the 
other two ensembles we conclude that gauge field config- 
urations made up solely of akyrons are not suited to pro- 
duce confinement. On the other hand, gauge field con- 
figurations, which are responsible for confinement, nec- 
essarily contain instantons and antiinstantons. 

It is interesting to note that any linear superposition 
of akyrons and, therefore, any field configuration in a 
pure akyron ensemble have zero topological charge den- 
sity (c.f. Appendix |B|) . Since akyron ensembles do not 
exhibit confinement, this supports the common expecta- 
tion that confinement and topological charge are closely 
related. 



V. SUMMARY AND OUTLOOK 

In this work we have presented the pseudoparticle ap- 
proach, a numerical method to compute path integrals 



in effective SU(2) Yang-Mills theories. We have calcu- 
lated the static quark antiquark potential at zero and 
at finite temperature, the topological susceptibility and 
the critical temperature of the confinement deconfine- 
ment phase transition in different pseudoparticle ensem- 
bles. The pseudoparticle approach is able to reproduce 
many essential features of SU(2) Yang-Mills theory with 
a comparatively small number of degrees of freedom. 



A. The pseudoparticle approach as a successful 
effective theory 

When using 400 instantons, antiinstantons and aky- 
rons, the static quark antiquark potential is linear 
for large separations with an attractive l/i?-correction, 
which is in qualitative agreement with the bosonic string 
picture and with results from lattice calculations. The 
string tension a is not only positive but also an increas- 
ing function of the coupling constant g. Therefore, when 
the scale is set by the string tension, one can adjust the 
size of the spacetime hypersphere in physical units by 
choosing an appropriate value for g. 

We have also calculated the the topological suscepti- 
bility x an d the critical temperature T cr jticai- The di- 
mensionless ratios y^l^ja 1 ^ 2 and Td-iticai/c 1 / 2 are con- 
stant for a wide range of coupling constants, i.e. a. \ 
and T cr ; tica i exhibit consistent scaling behaviors with re- 
spect to g. This success strongly indicates that the pseu- 
doparticle approach has the potential to reproduce cor- 
rect Yang-Mills physics. The values of both dimension- 
less ratios are of the right order of magnitude compared 
to lattice results. 



B. Properties of confining gauge field 
configurations 

For ensembles made up of instantons, antiinstantons 
and akyrons the string tension shows only a weak depen- 
dence on the pseudoparticle size A. It seems that confine- 
ment is mainly a consequence of the long range behavior 
of the building blocks, which is unaffected by A. This has 
been confirmed by considering ensembles of Gaussian lo- 
calized pscudoparticles and ensembles of singular gauge 
pseudoparticles, for which the size parameter A strongly 
affects the long range behavior. For small A there is 
only little overlap between neighboring pseudoparticles 
and there is no sign of confinement. Increasing A to a 
value, where pscudoparticles overlap and interact signif- 
icantly, restores quark confinement. The conclusion is 
that gauge field configurations, which arc responsible for 
confinement, contain extended structures. On the other 
hand, gauge field configurations with only localized exci- 
tations do not produce confinement. 

Comparing our "standard ensemble" with a pure in- 
stanton and antiinstanton ensemble and a pure akyron 
ensemble has shown that confinement arises due to in- 
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stantons and antiinstantons and not because of akyrons. 
Keeping in mind that gauge field configurations made 
up solely of akyrons have vanishing topological charge 
density, our findings support the common expectation 
that topological charge and confinement are closely re- 
lated. For quantitative results akyrons seem to play an 
important role. The dimensionless ratio x 1 / 4 /^ 1 / 2 is sig- 
nificantly closer to the lattice result, when there are not 
only instantons and antiinstantons but also akyrons. 



C. Outlook 

An obvious possibility for future research is to calculate 
correlation functions in order to extract glueball masses. 
Comparing the resulting masses with results from lattice 
calculations would be another check of the pseudopar- 
ticle approach. Furthermore we plan to compute other 
observables at finite temperature, e.g. the energy density 
and the pressure. 

A major new direction is to include fcrmions in the 
pseudopartiele approach. The goal would be to obtain a 
model for SU(2) Yang-Mills theory, which exhibits both 
chiral symmetry breaking and a confinement deconfine- 
ment phase transition. 
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APPENDIX B: ANY LINEAR SUPERPOSITION 
OF AKYRONS HAS VANISHING 
TOPOLOGICAL CHARGE DENSITY 

The gauge field of a single akyron with index i, ampli- 
tude A(i), color orientation matrix C ab (i) and position 
z(i) is given by 

= A(i)C ab (i)a b akymn Jx-z(z)) = 



A{i)C al (i) " 



{x-z{i)) 2 + X 2 



(Bl) 



(c.f. ©). In the corresponding field strength the deriva- 
tive terms cancel each other: 



= WB-MJB^fiKW. (B2) 



For any linear superposition of akyrons 



4" 



the same is true: 



For the topological charge density follows 



(B3) 



(B4) 
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APPENDIX A: COLOR ORIENTATION 
MATRICES 

A color orientation is a spacetime independent gauge 
transformation. It can be specified by an element of 
SU(2): U = c + ic a <j a with c 2 , + c 2 = 1. Applying 
such a gauge transformation to a gauge field yields 



K = ua.u- 1 

or expressed in components 
< = Tr(a a A') = 



(Al) A°(x) 



2c a c b + e abc 2c c c ) A 



(A2) 



=c* b 



We refer to C ab as color orientation matrix. Color orien- 
tation matrices fulfill 



CC T = C T C = 1 
det(C) = 1. 

Therefore, they are elements of SO (3). 



(A3) 
(A4) 



APPENDIX C: INSTANTONS AND 
ANTIINSTANTONS FORM TRANSVERSE 
GAUGE FIELDS, AKYRONS FORM 
LONGITUDINAL GAUGE FIELDS 

Any gauge field A® can be written as a sum of plane 
ives: 

1 



d 4 ke~* kx 3 



(CI) 



(2^) 4 

where A a ^, the Fourier transform of A^, is given by 
A;(k) = [ d 4 ke ikx A°(x). (C2) 



The Fourier transformed gauge field A c ^ can be decom- 
posed in a transverse and a longitudinal part: 



with 



A?. 



/^transverse 



(k) + ^longitudinal W (^3) 



A H, transverse (^0 
^.longitudinal^'') x ^M' 



(C4) 
(C5) 
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Superpositions of instantons J7J) and antiinstantons ((SJ) 
form transverse gauge fields, whereas superpositions of 
akyrons ((9|) form longitudinal gauge fields. This can be 
seen by considering the Fourier transforms of these pseu- 
doparticles: 



^i.instantonW — V^ v k v f(\k\) 
^antiinstantonW = Vpvkvf(\k\) 
a^akyron(fc) = ^Vd^l). 

where 



/(*) 



8ir 2 i f \k\XKii\k\X) /c 2 A 2 ^(|/c|A) 
W 



(C6) 
(C7) 

(C8) 
(C9) 



(K\ is a modified Bessel function of imaginary argu- 
ment). ([Co]) and (|C7|) satisfy (|C4|) due to the antisym- 
metry of 77" and 77^, while (|C8[) obviously fulfills (|C5[) . 



APPENDIX D: INSTANTONS AND AKYRONS 
FORM A BASIS OF ALL GAUGE FIELD 
CONFIGURATIONS 

In this appendix we show that "almost any gauge field 
configuration" can be represented by a linear superposi- 
tion of infinitely many instantons and akyrons. 

The starting point is the "continuum limit" of {T! 
without antiinstantons: 



j(.T — z) + 



z)C ab (j, z)< akyron - z) 

J=10 / 



(Dl) 



(the sum over all pseudoparticles has been replaced by 
an integration over spacetime; furthermore, nine instan- 
tons and three akyrons are allowed to share the same 
position). 

Inserting and © in (|Dl)l yields 



12 



J2A(j,z)C ab (j,z)S bl 5 lJtv 
i=io j 

=S»°(*) 



(a; - z) 2 + A 2 ' 



(D2) 



It can be shown that in general nine color orientation 
matrices form a basis of all 3 x 3-matrices Q. Therefore, 
any S ab and S a0 can be realized by suitably chosen am- 
plitudes A{i, z) and A(J, z). Hence, the problem has been 
reduced to the question whether any gauge filed configu- 
ration A" can be represented by suitably chosen S ab and 
S a0 . 



Fourier transforming (|D2[) turns the convolution into 
an ordinary multiplication: 



'S ab (k) v l+S a0 (k)5^ 



j4 

axe 



ikx 



X 1 + A 2 



= (S ab (k) v l„+S a0 (k)5^)k„f(k) 

(/ is defined by (IC9l) ). 
Without loss of generality we consider a = 1: 



A\{k) 
A\{k) 
\A\{k)J 



(D3) 
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fc -k\ -k 2 -k 3 

k\ fc -k 3 k 2 

ki k 3 k -ki 

h ~k 2 ki k 

V 

=/c(&) 



/5 io (fc)\ 

5 u (fc) 
S 12 {k) 
\S 13 (k)J 



(D4) 



For k ^ this equation can be solved for 
(S w ,S n ,S 12 ,S 13 ), because / ^ and 
det(/C) = |fc| 4 7^ 0. For fc = both / and K are singular. 
To study this case, we first deduce 



d 4 xA a Jx) = 



(D5) 



from (|D2j) by applying a proper regularization scheme. 
(|D5j) implies = 0) = 0. Inserting this in (fD4|) 

shows that the value of S aB (k = 0) has no effect on the 
gauge field A". On the other hand, changing S aB (k = 0) 
amounts to a constant shift of S aB : S aB -> <S aB + <S£ B . 
That is adding Sq B changes S aB , whereas the gauge field 
A® remains unaltered. To get rid of this redundancy, we 
require 



(DO) 



To be able to represent any gauge field configuration A^ , 
we have to find a way around (|D5[) . This can easily be 
achieved by adding constants to (|D1[) . 

The final result is the following: any gauge field con- 
figuration A a ^ has a unique expansion 

A a Ax) = 



J ^ z (^2A{i,z)C a \hz)a 



/j.instanton 



(x — z) + 



12 



J2 A(J, z)C ab (j, z)< akyron 0r -z)\+B; (D7) 
j=io / 

in terms of ,4(j,z) and -4(j, z) constrained by (|D6|) and 
B°, where C ah (i,z), i = 1, . . .,9, as well as C al (j,z), 
j = 10, . . . , 12 arc linearly independent color orientation 
matrices. 
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